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Abstract 

Given an efficient and systematic metliod for generating input sets for free fermionic 
heterotic model building we consider what the realistic bounds are for a statistical anal- 
ysis of the free fermionic Landscape with a classical computer. We then consider what 
kind of improvement could be expected on a quantum computer. We do outline the ba- 
sic structure of the relevant quantum algorithms, but we do not detail the construction 
of the oracle which would be involved in the calculations. 



* m_robinson@baylor . edu 
^gerald_cleaver@baylor.edu 



1 Free Fermionic Models and the String Landscape 

During the decade following the first string revolution, the focus of string phenomenol- 
ogy was two-fold: development of methods for constructing consistent string models with 
compactified dimensions and searching within the domain of each construction method for 
those models with at least quasi-realistic phenomenology pQ. Many three generation string 
models were found within a few years. The dominant view back then was that discovery of 
the true string vacuum (model) was within reach-that the string "needle" would eventually 
be found within the stringstack of but a few trillion vacua. Eventually a handful of quite 
realistic MSSM or Near-MMSM three generation models were indeed found |2], especially 
following the first Minimal Standard Heterotic String Model [3] . 

However, following the second string revolution the rise of M-theory has taught string 
phenomenologists the likely impossibility of finding a "true" string vacuum somewhere on 
the string/M-landscape composed of at least 10^"° vacua. In addition to the vast number 
of vacua, all vacua now appear to be on equal footing. Thus, the phenomenological goal 
has shifted from studying individual string models to better understanding statistically the 
characteristics of the string/M-models on the landscape, or at least within specific domains 

The (often overlapping) domains on the landscape frequently correspond to model con- 
struction methods. One construction method that has been widely explored in terms of 
individual models and for which large scale statistical studies are underway is free fermions 
O E]. The free fermionic heterotic string has provided many quasi-realistic (Near-)MSSM- 
like models semi-GUT models [TOllII], and GUT models |12j. In the context of the 

second string revolution, the important features to now investigate are the overall phenemo- 
logical patterns within these and the many more, as yet undiscovered, free fermionic models. 
Several studies of the statistical properties and characteristics of free fermionic heterotic 
models are underway by different research groups [131 HH UHl [151 [13 [IB] • 

An algorithm to systematically and efficiently generate free fermionic heterotic models 
was recently introduced [15]. As a first application, we have initiated an indepth study of 
the statistics of the gauge groups in free fermionic heterotic strings. Our approach enables 
a complete study of all gauge group models to be generated and analyzed with extreme 
efficiency, up to continually increasing Layers (the number of gauge basis vectors) and Or- 
ders (the lowest positive integer multiple N that transforms each basis vector back into the 
untwisted sector mod(2)). Our algorithm significantly increases by several orders of mag- 
nitude the efficiency of a systematic statistical study of regions of the parameter space of 
free fermionic heterotic models [15]. Nevertheless, current computational capabilities pose a 
serious limit to the range of parameter spaces that can be systematically investigated. 

Thus, given an efficient and systematic method for generating the required input data for 
a given string model construction method, such as the free fermionic heterotic approach, we 
consider what the realistic bounds are for a statistical analysis of the Landscape with a classi- 
cal computer. We then consider what kind of improvement could be expected on a quantum 
computer. We do outline the basic structure of the relevant quantum algorithms, but we 
do not detail the construction of the oracle which would be involved in the calculations. 
Since we will be relating quantum computing techniques to free fermionic construction, we 
briefly review in Section Two free fermionic heterotic model building and modular invariance 
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constraints. Then in Section Three we review the general concepts of quantum computing. 
In Section Four we discuss Grover's Quantum Search Algorithm, followed in Section Five by 
its application to the statistics of the free fermionic Landscape. Section Six concludes our 
discussions. 



2 Free Fermionic Heterotic Models 



The first object required to specify a model in the free fermionic heterotic string [5l E] is 
a set 

A = {a' e Q'^ n (-1, in e {1, L G N}}, (2.1) 

where components a], j = 1, . . . , 20 are boundary conditions for real worldsheet free fermion 
degrees of freedom of the left-moving supersymmetric string, and a* for j = 21, . . . , 64 are 
boundary conditions for real worldsheet free fermion degrees of freedom on the right moving 
bosonic string. In the C basis (complex fermions), each component of a* is double counted, 
and a* is a 32 (10 + 22) component vector (which can be generalized to include left-right 
paired real fermions). 

The order Ni of a given a* is defined agj 



Ni = mm{m G N | ma] = mod 2 Vj}, (2.2) 

with 

Nij = LCM{Ni,Nj). (2.3) 
Thus, each component of an a of order is of the form 

(a). G^n[-l,l) (2.4) 

Modular Invariance demands 

NiO^ -0^ = mod 8 (2.5) 

Ni- a' ■ a^' = mod 4 (2.6) 

in the C basis [3, |6] , along with requiring that the number of real rermions simultaneously 
periodic for any three basis vectors is even (including cases where one or more of the basis 
vectors may be the same). Further, each model contain = I, the 64 real-component vector 
with every element equal to one. 

In a given free fermionic model, the different sectors of a model are formed by all linear 
combinations of the a*'s, with coefficients G N where each coefficient < N^. Each 
linear combination, or sector, is denoted 

L 

V'' = J2mi mod 2. (2.7) 

i=l 



*Wc include "0" in the set of natural numbers N as in set theory, in contrast to its general exclusion in 
number theory. 
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We can (and will) think of each set of coefficients as an L- dimensional vector in N^, 
whose i*^^ component is constrained by the order of a*. 

For a given sector V'^, a worldsheet fermion transforms as 



/. ^ _e-^ /. (2.8) 

around non-contractible loops on the worldsheet. Thus, for M fermions, Vj' must be either 
or 1, whereas for C fermions, Vj' must be rational. For each sector, we can form the U{1) 
charges for the Cartan generators of the unbroken gauge groups (which are in one to one 
correspondence with the U{1) currents f*fj for each complex fermion /,); 

Q^, ^ ^V' + (2.9) 

where is a fermion number operator which counts each mode of fj once and of f* minus 
once. Or, in other words, F.'^g{ — l,0,l}Vi. 

The second object required to specify a model in the free fermionic heterotic string is an 
L X L matrix kij. Modular Invariance imposes the following constraints on kij; 

kij + kji = ^cP ■ mod 2 (2.10) 
1 • 

hi + kit = -a* ■ a* - Sj mod 2 (2.11) 

Njkij = Omod 2, (2.12) 

where Si is the 4 dimensional spacetime component of a*. Furthermore, kij G (—1, 1]. The 
GSO projection constraint for physical states is; 

L 

■Qv>==Yl ^nhn + Si mod 2. (2.13) 

n=l 

This completes our review of the general construction the free fermionic heterotic string 
models. We now now consider how quantum computing could significantly aid a systematic 
statistical study of this region (or of any other region) of the Landscape. 



3 Quantum Computation 

In the years since 1965 when Gordon Moore proposed it, Moore's Law has approximately 
held true. However, due to the drastic decreases in processor size and the increasing effects 
of quantum mechanics, it appears that fundamental fabrication barriers lie in the not too 
distant future]^ The physical construction of classical computers seems to demand a "smallest 
scale" beyond which we cannot go, and it is likely that we will reach this scale in the next 
decade or two. 

One solution to this problem is to reformulate the physical construction of computers. 
One such paradigm is provided by Quantum Computation. The basic idea of Quantum 

^The content of this introduction comes primarily from ^19j and [20] 
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Computation is that that the smallest unit of information is not a bit, but rather is a 
quantum bit, or q-bit, which can exist in a superposition of 1 and at the same time; 

|^) = a|0)+/3|l), 

where |ap + |/3p = l, and |1) = and are the computational basis vectors. 

There then exists a set of universal quantum gates which, similarly to the classical uni- 
versal gates, can be used to build an arbitrary operation on a set of q-bits. As an example of 
why such a system is useful, consider one such universal quantum gate, the Hadamard gate 

acting on each q-bit in an n q-bit register |0)'^"' = |0)i|0)2 ■ • • |0)„. This will give 

^ 2"-l 



where i runs over all 2" n q-bit states. With a classical computer, we could act on a single 
n bit register with some operator C, taking C\i)®^ — l^')*^", and this operator must be 
used once for every possible input. However, with the quantum computer, we can act on 



the superposition (3.2) a single time with some operator U (we use U because quantum 



operators must be Unitary), and 

^ 2"-l 

will be a superposition of every outcome at once. So using only n q-bits, and a single use of 
f/, we have 2" transformations. 

While this is impressive, there is a serious drawback. While this paradigm offers a massive 
parallelism, we obviously don't have access to every state in the superposition. All we can 
do is make a measurement of each q-bit, and we have equal probability of getting any n q-bit 
state. Getting a practical result out of a quantum computer therefore takes quite a bit more 
work. We mention this example merely to give an idea of how quantum mechanics could be 
used in computation. 

Despite decades of work, there has been relatively little progress made in developing 
algorithms for quantum computers. One problem is that a quantum algorithm is only ac- 
ceptable if it offers a speedup over the fastest classical algorithms. Another problem is that 
all of our intuitions come from the classical (looking) world we experience, and therefore 
we lack clarity in "thinking quantum" in order to design algorithms. However, there are a 
few quantum algorithms that have been developed and which offer a significant increase in 
efficiency over their classical counterparts. These existing algorithms could potentially be 
used in generating statistics on the free fermionic heterotic Landscape. We will therefore 
outline how these algorithms work, and discuss what type of improvement they would offer 
in such an analysis. To date, there are no quantum computers in active operation. However, 
it is likely that in the next few decades their production will allow for applications like the 
one discussed in this paper. 
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4 Quantum Algorithms 

For a more detailed exposition of the nature of quantum gates and their physical con- 
struction, see [I9], which we follow closely in this section. It should be assumed that all of 
the gates discussed below are realizable. 

4.1 Quantum Fourier Transform 

The first quantum operation we consider is the quantum Fourier transform (QFT) [TU]. 
The action of the QFT on an n q-bit state is defined by 



2"-l 

But we can re-express this as 

2"-l 



^ 2"-l 



br - 4^^^"^!^)^" (4-2 

(|0) + e2^^o->^|l)) (|0) + e2^^o.j„_ij„^ . . . (||o) + e2^io.jii2-j„ 



^2" 



(4.3) 



where O.jkjk+i ■ ■ ■ jm = ^ + + ■ ■ ■ 2m''-^+i is the binary fraction for j, where each ji is 
defined by j = ji2"-i + j22""2 + ■ ■ ■ + jn2^. 

4.2 Phase Estimation 

Suppose we have some unitary operator U with eigenvector \u) and eigenvector e'^'^^'^, 
where (p is unknown. We can use the QFT (4.3) to estimate to arbitrary accuracy [19j. 
We assume that we have oracles which will both prepare the state \u) as well as perform 
a controlled- 1/^^ operation for non-negative integers j. We then prepare a register with n 
q-bits all in the state |0) (10)®""), where n depends on how accurately we want to estimate 
0, and a second register in the state \u), or 



|0)®"|m). (4.4) 



We act on each of the q-bits in the first register with a Hadamard gate (3.1), creating 
the state 

+ |1))(|0) + |1))---(|0) + |1)) \u). (4.5) 



Next, we act on the second register, \u), with a controlled-f/ operation with U raised to 
successive powers of 2, creating the state 

^(|0) + e2-2"-V|i)) (|o) + e^-^-^-^li)) . . . (|o) + e2-2V|i)) |^) 



2"-l 

_^e(^-)/c<^|fc)H. (4.6) 



k=0 



Finally, we can do an inverse QFT to the first register, taking 



2"-l 

e(2-)^^|fc)|«)^|0)|«), (4.7) 



'On .^^^ 



^ fc=0 



where is an estimator for when the first register is measured. If n q-bits are used in the 
first register (as indicated above), then this algorithm will accurately approximate to 



n 



log2 ( 2 + ^ 



(4. 



bits with probability 1 — e. 

4.3 Grover's Quantum Search Algorithm 

We now consider the algorithm at the heart of our application to superstring model 
building statistics, Grover's search algorithm [19] [20] [21]. The idea is that we start with 
a space of = 2"^ elements which we want to search through. We index the elements in 
the space, which requires n bits. We are looking for elements which satisfy some specific 
criterion, of which assume there are M in the search space. 

Next, we define what is called an oracle, which is a function f{x) which takes an n-bit 
value X G {0, ... — 1}, and has the value f{x) = if x is not one of the M solutions, and 
the value f{x) = liix is one of the M solutions. We will not concern ourselves with the 
construction of the oracle in this section or the next. We simply assume its existence. 

The actual operation of the oracle is as follows; given an n + 1 q-bit state prepared as 
|a;)®"|?/), the oracle operator O makes the transformation 

0\x)®''\Q) = \x)®''\y® f{x)) (4.9) 

(where © denotes addition mod 2). In other words, O changes the value of the last q-bit iff 
X is a solution. 

A more useful application of this is to prepare the last q-bit in the state |0) and do a 



Hadamard transformation (3.1). Defining |— ) = :^(|0) — |1)), we then have 

C)|x)^"|-) = (-l)^(")|x)|-). (4.10) 

To carry out the actual search algorithm, the computer is prepared in the state |0)®"|— ), 
and the each of the first n states are acted on by Hadamard gates, producing 



2" — 1 



x=0 

To this state we apply the 'Grover Operator" G defined as 

G = if®"(2|0)(0| - /)iJ®"0 = (2|V)(V'| - I)C- (4.12) 
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Let I a) be the normalized superposition of all states which are not solutions, |a) = 
T^lf Exs/(aO=ok)' and 1/3) those that are, \(3) = ^ Exs/(aO=i 1^)- Then, ( |4lT| ) can 
be written as 




(we have dropped |— ) for notational simplicity). 

It is clear that C is a reflection in the (|a), |/3)) plane about \a), and that 2\ip){ip\ — I is 
a reflection in the same plane about So, G is a rotation in the (|a), |/3)) plane. If we 
define 



2 \ N 

then k successive applications of G will produce the state 



cos- = \/ — - — , (4.14) 



G"^|^) = cos ( ^^T^^l |a) + sin f ^^^^^ ) (4.15) 



If we do approximately k = \\J~^ applications of G, we will have rotated the state very 
near so that measurement of will give a solution with high probability. So, whereas 
with a classical computer, 0{N/M) applications of the oracle are needed, only 0{^^ N / M) 
are needed with the quantum computer. 

4.4 Quantum Counting 

Finally we review the Quantum Counting algorithm [19] , which combines ideas from the 
previous sections. Grover's algorithm required that we know both the total search space iV, 
as well as the number of solutions M . But, it is possible that M is not known in advance. In 
this case, there is a fast quantum algorithm which will calculate M to high probability. The 



primary idea is to use the phase estimation algorithm of section 4.2 to find the eigenvalue 
of the Grover operator G, from which we can calculate 9 (cf. section 4.3), and therefore M 
(assuming we know N). The eigenvalues of the Grover operator G in the space spanned by 
|a) and \(3) are e^*^. We can therefore use the phase estimation algorithm to find an estimate 



for 9 to whatever accuracy we wish. Once 9 is obtained, we can use (4.14) to find M. 

5 Statistics on the Free Fermionic Landscape 
In [T3] we found that for a given layer L and total order product 

L 

A^\[{m-i, (5.1) 

i=l 

where Ni are the respective orders of the i = 1, . . . L a*'s, the number of modular invariant 
input sets can be approximated by 

L i-1 
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Collecting statistics then consists of analyzing each model, looking for occurrences of a 
particular gauge group. In other words, if we want to know what percentage of all layer 
L and order A sets contain a group M, we must analyze all B sets. So, given an oracle 
which takes a set of a*'s and returns 1 if M is produced and otherwise, the oracle must 
be invoked B times to get an exact statistic. Obviously there will be a limit on A and L, 
beyond which we can no longer perform comprehensive investigations. Graphing (5.2) vs. A 
logarithmically (at L = 1), 
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we see that the computational limit will be between 20 and 25. 



With a quantum computer, however, we can combine Grover's algorithm (section 4.3) and 



the phase estimation algorithm (section [4.2 ) to improve our search capabilities drastically. 



Equation (5.2) gives us a good approximation for the total number of possibilities at a given 



L and A, and therefore all we need to know is the phase of the Grover operator eigenvalue 
for an oracle that returns 1 for group M, and we can calculate the percentage of sets (at 
that L and A) that will contain M . 



For L and A, we can index each set using at most log2(-B) q-bits. So, equation (4.8) 
demands that the first register in our phase estimation algorithm contain 



k = log2 B + 



loga 



1 

Ye 



(5.3) 



and therefore the phase estimation will involve only k invocations of the oracle. So, if we 
want to measure the eigenvalue with 99% probability, our phase estimation procedure will 
only involve 



k = lo£ 



L i-l 

i=i j=i 



\ r(23 + A) 



7 2M!(A + l)^r(23) 
We can graph k vs. A logarithmically (at L = 1), 



+ 6. 



(5.4) 
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5 10 J 5 20 25 30 35 40 

which is obviously a drastic improvement over the classical approach. We can see the im- 
provement more clearly by graphing js- 




\Q 15 2fl 25 30 35 40 



6 Conclusion 

When Grovers search algorithm is combined with the phase finding algorithm based on 
the quantum Fourier transform, an enormous improvement in the depth of statistical analysis 
on the free fermionic heterotic string Landscape over classical computation can be obtained. 
The massive parallelism made possible by the superposition of states would render such a 
statistical search effectively unbounded. However, because quantum computing is still in its 
infancy and no quantum computers are in operation, the purpose of this paper is merely 
to outline a specific application of an already existing quantum computing algorithm for a 
potential future use in string theory. 
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